clc;
clear;


tspan=[0,40];	%求解的时间区间
y0=[1;0];			%初值向量应与DyDt.m文件中y形式一致。
[tt,yy]=ode45(@DyDt,tspan,y0);					%<3>
plot(tt,yy(:,1))
xlabel('t'),title('x(t)')  
